% Amplitude of current for odn plazma
% Usage: y = sum_I_odn_pl(phi, CONSTS)

function y = sum_I_odn_pl(phi, CONSTS)
    k0     = CONSTS.k0;
    a      = CONSTS.a;
    c      = CONSTS.c;
    d      = CONSTS.d;
    phi_0  = CONSTS.phi_0;
    vacuum = CONSTS.vacuum;
    
    if (vacuum)
        h = k0;
    else
        h = CONSTS.h;
    end
    
    a0 = - 1i*c/(4*pi*k0*a*(log(4*a/d)));
    y = a0 * pi*h*a * ((cos((pi - abs(phi - phi_0))*h*a)) ./ (sin(pi*h*a)));    
end